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ABSTRACT 

Cosmic strings with small-scale structure have a coarse-grained mass per unit 
length jj, which is larger than the string tension. This leads to an effective Newto- 
nian gravitational line source and to a characteristic translational velocity which 
is smaller than for strings without small-scale structure. Here, the accretion of hot 
dark matter onto such strings is studied by means of the Zel'dovich approxima- 
tion. We find that clustering is greatly enhanced by the Newtonian line source. In 
the limit of vanishing translational velocity, the first nonlinear filaments form at a 
redshift of greater than 100 for standard values of \x. 



1. Introduction 



A theory in which the dark matter is hot and in which the primordial density 
perturbations are seeded by cosmic strings is a viable alternative 1 ' 2 -* to the currently 
popular variants of cold dark matter (CDM) cosmologies with adiabatic density 
fluctuations produced during a hypothetical period of inflation in the very early 
Universe. On large length scales, the theory predicts a scale-invariant spectrum of 
perturbations. The normalizations of the model from large-scale structure 3 ) and 
from cosmic microwave background anisotropy measurements 4 -* are in encouraging 
agreement. The theory also predicts that on scales comparable to and larger than 
the comoving horizon at t eq , the time of equal matter and radiation, galaxies are 
concentrated in sheets or filaments 3,5-8 ) , in agreement with the results from the 
Center for Astrophysics galaxy redshift survey 9 ) . 

If the primordial density inhomogeneities are nonadiabatic, slowly decaying 
seeds such as cosmic strings, hot dark matter (HDM) is a viable dark matter 
candidate. In contrast to adiabatic inflationary models in which free streaming 
erases perturbations on all scales smaller than the maximal free streaming length 
X l j ax , the strings survive neutrino free streaming 10 ) and can at late times seed 
structures on scales much smaller than Xj ax . There is 11 ) still power on scales 
smaller than \™ ax (although less than in the inflation-based CDM cosmology), and 
the power spectrum may well be in better agreement with the results from recent 
observations 12 ) (the decrease in power on small scales compared to the inflationary 
CDM model is a promising feature). 

In spite of the abovementioned advantages of the cosmic string theory, there is 
still a lot of ignorance about the specific predictions of the model. One important 
reason is that the initial conditions for structure formation in the theory are not 
yet accurately known. It is well established that the network of cosmic strings 
(and thus the distribution of the seeds for density inhomogeneities) approaches 
a scale-invariant distribution, but the detailed properties of this distribution are 
uncertain. It can be shown that at all times t, there are a few strings which are 
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straight on a correlation length £(£) ~ t crossing each Hubble volume, but the 
small-scale structure of these strings is not known. Also uncertain is the exact 
number of such strings per Hubble volume. 

The evolution of the cosmic string network may leave behind on the strings a 
lot of structure (kinks and wiggles) on scales smaller than £. The typical ampli- 
tude of these perturbations is much smaller than £, and hence in a coarse grained 
description, the strings can still be viewed as straight on a scale £, but with renor- 
malized mass per unit length and tension. Since the integrated length of the string 
between two points is longer than the separation between the points, the coarse 
grained mass per unit length is larger than the microscopic length. At the same 
time, the coarse grained tension is less than the microscopic tension. As first 
discussed by Carter 13 ) and subsequently analyzed in detail by Vollick 6 ) and by 
Vachaspati and Vilenkin 7 ) , small-scale structure on strings can lead to substantial 
modifications of the motion and gravitational effects of the strings. 

A long straight string without small-scale structure induces no local gravita- 
tional force. However, space perpendicular to the string is a cone with deficit 
angle 14 ) 

a = 8vrG/i, (1.1) 

where \i is the mass per unit length of the string coarse grained on a scale £ and 
G is Newton's constant. When moving with a transverse speed v s , such a string 
will induce a velocity perturbation of magnitude 

Sv = AnGpLVs^fs (1.2) 

towards the plane behind the string, leading to the formation of a "wake" , a planar 
density enhancement 5 -* . In the above, 7 S is the relativistic 7 factor associated with 

Vs- 

A long straight string with small-scale structure will induce, in addition to 
the deficit angle given by Eq. (1.1), a nonvanishing local gravitational force of 
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magnitude 6 ' 7 ^ 



F = 2mG(fx-T)/r 



(1.3) 



towards the string when acting on a test particle of mass m located a distance r 
from the string. The string tension is T, and T is related to /i by 13 ) 

liT = M §, (1.4) 

where /io is the microscopic mass per unit length of the string. In addition, when 
coarse grained on a scale of £, strings with /i > /iq have a smaller transverse 
velocity than those with \i = /iq. In the limit when the distance the string moves 
in one Hubble expansion time is small compared to the thickness of the wake, the 
structure seeded by the string becomes filamentary. 

It is important to know at what redshift the first density perturbations become 
nonlinear (this gives an upper bound to the redshift of first star formation). To 
answer this question in our theory, the clustering of HDM induced by cosmic strings 
must be studied. 

The accretion of HDM onto cosmic string loops was studied in detail in Ref. 
10. It was shown that free streaming slows down but never reverses the growth of 
perturbations on small scales. Structures form "inside-out" in the sense that the 
first shells to go nonlinear around the seed are the innermost ones. In a toy model 
in which only cosmic string loops are present as seeds, this leads to a "bottom-up" 
scenario for structure formation. The first nonlinear objects form at high redshifts. 

However, the current cosmic string simulations 15 ) indicate that most of the 
mass in the string network resides in the long strings (strings with curvature radius 
£(i) ~ t at cosmic time t). Hence, in order to determine the earliest time at which 
large-scale nonlinear structures can form, it is necessary to study the clustering of 
HDM induced by long strings. This was done in Ref. 3 for strings without small- 
scale structure. The results were very different than for loops. Planar collapse of 
HDM does not proceed in an "inside-out" fashion. Rather, the first sheet to go 
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nonlinear has an initial comoving distance q ma x ~ 1 Mpc from the center of the 
wake, and the onset of nonlinearity occurs late, i.e. at a redshift z max ~ 1 for 
values G/i ~ 1CT 6 (these results will be reviewed at the end of Section 2). 

In this paper, we study the accretion of HDM onto long strings with nonvanish- 
ing small-scale structure (which leads to the presence of a Newtonian gravitational 
line source). The calculations will be done in the limit in which the tranlational 
motion of the string can be neglected. 

We find that in this situation, structures once again form "inside-out" , and that 
the first nonlinear structures emerge at a high redshift (z max > 10 2 for G/i ~ 1(T 6 ). 
This result indicates that in a cosmic string theory with HDM, there will be no 
problem in explaining the recent observations of high redshift clusters 16 -* , quasars 17 ) 
and QSO absorption line systems 18 ) . 

In the following section we summarize our calculations, and in Section 3 we 
discuss some implications of our results. Units in which c = ks = 1 are used 
throughout. The scale factor of the Universe is denoted by a(t), t being cosmic 
time, and the associated redshift is z(t). The present time is to- 

2. Accretion of Hot Dark Matter onto a Line Source 

We will study the accretion of collisionless hot dark matter (HDM) onto a 
stationary cosmic string using an adaptation of the Zel'dovich approximation 19 ) to 
HDM. This modification was applied to planar collapse in Ref. 3 and was shown to 
give good agreement with the results of a linearized Boltzmann equation method 
(which takes into account the full phase space distribution of HDM particles; see 
Refs. 3 & 10). 

HDM particles have near-relativistic velocities v(t) at times close to t eq . This 
prevents them from clustering on length scales smaller than their free streaming 
distance 

Xj(t) ~ 3v(t)z(t)t : (2.1) 
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the average comoving distance a particle will move in a Hubble expansion time. 
Note that since v(t) oc z{t) oc a(£) -1 , the free streaming length decreases in time 
as 

Xj(t) oc t" 1/3 (2.2) 
and takes on its maximal value \™ ax at t eq . In a spatially flat HDM Universe 

v(t eq ) =v eq ~0.l (2.3) 

and hence 

ynax _ ^h^Mpc, (2.4) 

substantially larger than the typical size of a galaxy. 

We will first describe the Zel'dovich approximation for filamentary accretion 
and then in a second step adapt the method to HDM. Without loss of generality 
we can choose coordinates in which the string lies along the z coordinate axis. 
We consider the response of the dark matter particles to the gravitational line 
source produced by the small-scale structure of the string, focusing on the "plane" 
(strictly speaking a cone) at z — 0. Let q denote the initial comoving separation 
of a dark matter particle from the string. In the absence of the string line source, 
the physical separation f(t) would grow as the Universe expands. The presence of 
the string, however, leads to a comoving displacement ip(q,t) which grows in time: 

r(q,t)=a(t)(q-j(q,t)). (2.5) 

The Zel'dovich approximation consists of taking the Newtonian equations 

f=-JU(f,*) (2.6) 
for f(q,t) with the gravitational potential $(r, t) determined by the Poisson equa- 
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tion 



S2 * = 4n G ( P (r,t) + XS <- r > 



dr 1 



2nr 



(2.7) 



and linearizing in ip. On the right hand side of the Poisson equation, the total 
density has been separated into the dark matter density p(f, t) and the contribution 
of the string-induced line source with an effective mass per unit length 



A = /i-T 



(2.8) 



(see Eq. (1.3)), and where S(r) denotes a Dirac delta function in cylindrical co- 
ordinates. The density p(f, t) is in turn determined in terms of the background 
density po(t) by the mass conservation equation 



p(f,t)d 3 r = a 3 (t)p d 3 q 



(2.9) 



which gives 



p(f,t)=a 3 (t)podet 1 



dr 
dq 



a^po^r 1 ^ po(i + §^A(q,t)), (2.10) 



where r, q and ip r are (in cylindrical coordinates) the radial components of r, q 

—* — * 

and ip, respectively, and where in the final step we have linearized in ip. 
Inserting (2.10) into (2.7) and integrating once gives 



d 



dr 



-$ = 47TG 



1 - X9(r)r 

- PQ (t)f+po(t)a(t)iP(q,t) + 



(2.11) 



which we insert on the right hand side of (2.6). Note that f is the unit radial 
vector in the plane perdendicular to the string. On the left hand side of (2.6) we 
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substitute the second time derivative of (2.5), 

— r = aip + 2aip — dq + dip. (2-12) 



Using the FRW equation d = —^j^-poa (valid in the matter dominated era) to 
substitute for poa in (2.11), we obtain as our basic equation 

•■ n d ■ d , 2AG 

Vv + 2-A + S-iJr = — . 2.13 

This equation is only valid for t > t eq , since we made use of an approximate 
pressureless equation of state. Hence, we should also substitute a(t) = (t/to) 2 / 3 
and then obtain 

4 2 2XGt^^ 

Vv + f^r - 3*-=Vv = ?/ ^ y (2.14) 

as equation of motion for the radial displacement tp r . 

A second reason why the above derivation is only valid for t > t eq is that the 
presence of the uniform radiation bath has been neglected. This radiation would 
have to be considered in the Poisson equation (2.7) and would stall the growth of 
perturbations before t eq . Hence, for a cosmic string theory with cold dark matter, 
the initial conditions for solving (2.14) are 



^(ts)=^(ts) = 0, (2.15) 



with t s = t eq (for strings present at t eq ) and t s = t{ for strings "appearing" at a 
later time. 

For a static string network, it does not make sense to consider tj > t eq . How- 
ever, for a dynamical string network, U > t eq would apply to string segments 
whose curvature radius £ is equal to the Hubble radius at time U and thus start 
their proper motion at that time. 
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To first order in perturbation theory, we solve (2.14) using the Born approxi- 
mation 

v 4 _i ; 2 _ 2 , 2A4 /3 G 

A + - 3 t~ A - -t- a = —^r- 

This can be solved using the Green's function method which gives 



(2.16) 



t t 
Mt) = Mt) J dt'e(t f )u 2 (t')f(t f ) - u 2 (t) J dt f e(t') Ul (t f )f(t'), (2.17) 



where ui(t) = (|^) 2//3 and U2(t) = (^) 1 are the two fundamental solutions of the 
homogeneous equation, e(t) is the Wronskian 



e(t) = (u\U2 - U2U1) \ 



(2.18) 



and 



/(*) = 



2AG 



4/3 



Some manipulations give the end result: 



(2.19) 



ln(t/ts) - - I 1 - 



5/3 x 



[2.20) 



The real objective of the calculation, however, is to find the comoving scale q n \ 
which has gone non-linear at a given time t, as a function of the time of formation 
t{. This scale is defined as the scale at which at a time t the physical velocity of a 
particle r(q n i,t) vanishes, giving an equation for q n \ of 



a(t)(Qni(t s ,t)-> t p(t))=a(t)> t p(t). 



(2.21) 



Together with Eq. (2.20) and the assumption that a(t) = (t/to) 2 / 3 , this gives (after 
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a fair amount of algebra) 



9„i(*«.«)=« 1/3 «n /3 



i 



12GX 



5/3 



(2.22) 



For CDM, this equation gives the size of the non-linear region at time t for a 
structure which starts to grow at time t s . 

Because of neutrino free streaming, the above formalism is not directly appli- 
cable to the clustering of HDM. A rigorous analysis must follow the entire phase 
space distribution of the dark matter. Starting point is the collisionless Boltz- 
mann equation describing the time evolution of the phase space density, expanded 
perturbatively about a homogeneous thermal distribution. After integrating over 
momenta one obtains the Gilbert equation 20 ) , an integral equation for the con- 
figuration space density perturbation. However, in Ref. 3 it was shown that the 
Zel'dovich approximation can easily be adapted to HDM and then gives good 
agreement with the results from the Gilbert equation approach. 

The modified Zel'dovich approximation consists of using the same dynamical 
equation (2.14), but with initial conditions 



i/>{t 8 (q))=i/>(t 8 (q))=0, 



(2.23) 



where t s (q) is the time when the neutrino free streaming length Xj(t) equals q, i.e. 



*j(t a (q)) = q- 



(2.24) 



Solving for t s (q) gives 



ts — ■>■ 



Aq — 3v e qZ e q t e q. 



(2.25) 



Note that if (2.24) gives a value of t smaller than U, then t s (q) must again be 
chosen to be t{. 
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For HDM, turnaround in determined by Equation (2.22), but with the time t s 
depending on q as given by (2.25). Inserting (2.25) into (2.22) and solving for q at 
late times gives 



5z(t) 



( 



= 0. 



(2.26) 



This implicit equation gives, for a given scale q, the redshift z(t) at which this scale 
q first goes non-linear. This equation only holds for values of q smaller than qmaxi 
the value of q for which t s (q) equals ij, i.e. 



For q larger than q ma x, Equation (2.25) with t s = U must be used. 

Equation (2.26) can be solved numerically. The result is shown in Figure 1. 
The main conclusion is that in the filament model with HDM, non-linear regions 
can develop at very early times - at redshifts larger than 10 2 for the usual cosmic 
string normalization G/i = 1CT 6 . The first scales to go non-linear are small scales. 

These results contrast with those for planar collapse of HDM 3 ), where no non- 



linearities appear until a redshift of about 1, and the first scale to turn around is 
q = Qmax (see Figure 2). 



Qmax — 



(2.27) 
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q Figure 1: Redshift 

z n \ at which a scale q (in Mpc) turns around for filamentary accretion in the case 
of HDM (dashed line) and CDM (solid line). The value of GX was 5 • 1(T 7 , and 
the case ti — t e q was considered. 




q Figure 2: Redshift 

z n i versus scale q (in Mpc) for accretion of HDM onto a cosmic string wake for a 
value Gp, = KT 6 (from Ref. 3). 
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3. Discussion 



We have studied the accretion of HDM onto a stationary cosmic string with 
nonnegligible small-scale structure. Such a string acts as a Newtonian line source 
of gravity for HDM. Our main results are summarized in Figure 1: the first scales 
to turn around are small scales, and the redshift when the first non-linearities 
appear exceeds 10 2 for the usual cosmic string normalization. This contrasts with 
the results of planar collapse of HDM onto cosmic string-induced wakes, in which 
case no non-linear structures form until a reshift of about 1, and where large scales 
<? = Qmax are the first to turn around. 

These results are important in determining the redshift when the first non- 
linear filamentary structures form in a cosmic string model with HDM. The mass 
of a filament seeded at t eq by a horizon-length string is 

M(z) ~ p t eq z(t eq )q 2 nl (z) ~ 10 12 M o (M^) 2 , (3.1) 

where po is the background energy density and Mq is a solar mass. Such a mass 
is large enough to be the host of a quasar at a redshift much greater than l 21 ). 
Hence, explaining the origin of high redshift quasars and galaxies seems not to be 
a problem in this model. 

Obviously, much more work is needed in order to be able to reliably estimate 
the high redshift mass function in the cosmic string model with HDM. Besides the 
filamentary accretion mechanism discussed in this paper, there is also accretion 
at high redshift onto cosmic string loops 10 ). Which of the mechanisms is domi- 
nant unfortunately depends very sensitively on the parameters in the cosmic string 
scaling solutions which are not yet reliably known. 

Considering a long cosmic string as stationary in order to discuss the accretion 
of HDM on a given comoving scale q will only be a reasonable approximation if q is 
larger than the comoving distance the string moves at t{ in one Hubble expansion 
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time, i.e. 



q > t z(ti) 1/2 v s , 



(3.2) 



where v s is the transverse velocity of the string. For U = t eq , this gives 

q > 30v s Mpc. (3.3) 

For strings with v s oc 0.03, the approximation of filamentary accretion is hence 
reasonable for scales q oc IMpc. Note that the r.m.s. value of v s is about 0.15 in 
the current cosmic string simulations 15 ) which probably underestimate the effects 
of small-scale structure and thus overestimate v s . 

In conclusion, we believe our results are good news for the cosmic string model 
with HDM and should encourage further research. The first non-linear structures 
form early enough in order to explain the origin of high redshift quasars, galaxies 
and clusters. The results also imply that the value of G/j, could be lowered from 
its canonical value of Gfx ~ 10~ 6 without jeopardizing early structure formation. 
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